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ABSTRACT 


The tire plays a key part in many analyses used for design of aircraft landing 
gear. Examples include structural design of wheels, landing gear shimmy, brake 
whirl, chatter and squeal, complex combination of chatter and shimmy on main 
landing gear (MLG) systems, anti-skid performance, gear walk, and rough terrain 
loads and performance. This paper discusses which tire parameters are needed 
in the various analyses. 

Two tire models are discussed for shimmy analysis, the modified Moreland 
approach (ref. 1) and the Von Schl ippe-Dietrich approach (ref. 2). It is shown 
that the Moreland model can be derived from the Von Schlippe-Dietrich model by 
certain approximations. 

The remaining analysis areas are discussed in general terms and the tire 
parameters needed for each are identified. The conclusion of the paper is that 
accurate tire data allows more accurate design analysis and the correct pre- 
diction of dynamic performance of aircraft landing gear. 
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INTRODUCTION 


A number of important design analysis areas of aircraft landing gear require 
accurate data on tire dynamic parameters. One of the more obvious examples is wheel 
design. The primary loads acting on the wheel act through the tire onto the wheel 
bead seat area. The wheel cannot be efficiently designed without accurate data on 
tire to wheel contact loads for the entire spectrum of operating conditions to which 
the landing gear system is subjected. This includes landing impact, braked roll, 
taxiing, turning, cross wind landing, and in some cases obstacle impact. The latter 
is particularly true for carrier based airplanes. The wheels of such aircraft must 
accommodate running over arresting cables when the tire is in the near bottomed con- 
dition. Modern aircraft wheel design makes use of experimentally and analytically 
determined bead contact forces to arrive at optimized designs of the wheel cross 
sections. In spite of increased demands for more durable wheels having longer 
fatigue life and increased capabilities, wheel weight per unit load has decreased in 
recent years. 

Other design analysis areas involving the entire landing gear system include: 

• Landing gear shimmy (both nose and main landing gear) 

• Brake-landing gear whirl stability 

• Brake chatter and squeal 

• Complex combinations of chatter and shimmy 

• Anti -skid performance 

• Gear v/al k 

• Rough terrain operations 

This paper will concentrate primarily on the tire parameters needed in the first four 
of the above problem areas. 
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LANDING GEAR SHIMMY 


There are two tire models commonly used for aircraft landing gear shimmy 
analysis: the modified Moreland model (refs. 1, 3, and 4) and the Von Schlippe- 
Dietrich model (refs. 2, 5, and 6). The equations for the forces and moments at 
the tire footprint are similar for the currently used versions of the two models and 
as a result there are a large number of parameters which are common to both models. 
Both models require that the vertical load deflection characteristics of the tire be 
known along with the general load, inflation pressure and speed envelope of opera- 
tion of the aircraft. In addition the following parameters must be known for the 
operational envelope: 

m^ = mass of tire (kg) 

2 

Ipt ~ polar moment of inertia of inflated tire (kg*m ) 

2 

I^t “ diametral moment of inertia of inflated tire (kg-m ) 

= lateral spring rate of tire (N/m) 

C|^ = lateral equivalent viscous damping coefficient (N-sec/m) 

l-L^ = rolling tire torsional spring rate (N-m/m) 

I_Iq = rolling tire torsional damping coefficient (N-m/sec) 

Rg = loaded radius of the tire (m) 

= slope of drag force versus slip ratio for small slip ratios (N) 

R^ = tire rolling radius (m) 

For the modified Moreland model, two additional parameters are required: 

C = coefficient of yaw (rad/N) 

= yaw time constant (sec) 

These parameters are used in Moreland's cornering relationship: 


CF 


nj 


= ^ 


tj 






( 1 ) 
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where 


j = index indicating the jth wheel of a multi-wheel landing gear 
F . = lateral force acting on tire footprint (N) 

1 1 J 

= yaw angle or sideslip angle of tire footprint with respect to 
^ the tire (rad) 

Dots over a variable are used to indicate differentiation with respect to time. 

In the Von Schl ippe-Dietrich model it is shown that the cornering characteristics 
come about due to continuity between the tire distortion in the footprint's 
most forward position and the lateral distortion of the tire as a function of 
peripheral distance, combined with adhesion of the tire footprint centerline to the 
ground surface. Denoting the lateral displacement of the tire outside the foot- 
print as X and the peripheral distance around the tire as s (s = 0 corresponds to 
the forward tip of the footprint and s > 0 corresponds to moving upward around the 
tire center) it is found from experiments that when a pure lateral load is applied 
to the tire X is given by 


X = 



( 2 ) 


where Xp = lateral deflection of footprint (constant within footprint) (m) 
s = peripheral distance (m) 

a^ = experimental constant between zero and 1 (dimensionless) 

L = experimentally determined "relaxation" length (m) 

For type VII tires, aj is very close to 1. However, for many new design tires 
with a relatively small aspect ratio and a relatively small width between wheel 
flanges, a^ is significantly less than unity. For example, on a H45 x 17-20, PR22 
tire, typical values of a^ range from 0.80 to 0.92, depending on the vertical load 
and the tire pressure. 

As a consequence of Equation (2), it can be shown that when the wheel moves forward 
an amount ds the most forward position of the tire footprint will move an amount 


dX 


1 




(3) 


74 



I 


relative to its position prior to rolling forward. The quantity a^^/L is denoted by 

c = 3,-^/Jj (4) 

It can be shown that as a consequence of Equation (2) the steady state "cornering 
power" or cornering coefficient is given by 

N = % ^ (l + j (5) 

where 

N = cornering coefficient (M/rad) 

= tire coefficient defined by Equation (4), (1/m) 

h = half length of the tire footprint (m) 

For transient conditions the general kinematic situation is illustrated in Figure 1. 
Due to strut deflection, the wheel is tilted at an angle and twisted in the yavif 
direction at an angle xj/'s , The coordinate of the wheel axle center is ya and the 
coordinate of the most forward position of the tire footprint is y^p* The wheel is 
rolled forward an amount ds during which time the position variables of the tire 
change by amounts d(bL, d i// ' dy,, dX , and C^Xds as defined by Equations (3) and 

/ A \ D b ^ 



Figure 1.- Tire distortion nomenclature. 
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With this in mind we can, from Figure 1, write the quantity + dy^„ in two 
ways. We first consider the path a to b to c to d to e which gives^ ^ 


7fp + dy^p = Ya + Rg<^B + dy^ + Rgd0' - h(;^; + - (X + dX) (6) 


Next, we consider the path a to b to f to g to i to j which gives 


+ Rg0B - h 03 ' - X - + C^Xds (7) 


Equating (6) to (7) gives 

dy^ + Rgd0g - hdl//g'- dx = - %'ds + C^Xds 


( 8 ) 


which can be written as 

dX 


^ + C.X = T-^ + Rg — ^ + 0s h ^ 

ds X ds S ds ® 


(9) 


We 
from 


can eliminate X from Equation (9) by defining X in terms of y^ and y^ . Again 
m Figure 1, 

X = y^ - yfp + Rg0' - h0/ 


Substituting (10) into (9) gives 


dyfp 

~dT +^xyfp 


Cx(ya + Rg0s) - <1 + Cxh)0s' 


Or, letting -~ = D, this substitution will yield 
ds 


fp D +c 




- ( 


1 + C^h 


)^s'l 


( 11 ) 


( 12 ) 


As a result of adhesion of the tire footprint to the ground, we can write 



Taking the series expansion for the left hand side of Equation (13) gives 


yfp(s = -h) = yfp - h 


ds 


,2 


l£p 

2 1 ds^ 


(14) 


1 

Applying (14) to (12) and making the substitution D = — where = d/dt gives 


/ n »» 0 


3!v 


3-0^3 + . . .)(cx(y, + Vb> - + “ 

+'=x)k + ><»<#> B -^) 


(15) 


Equation (15) is the basic cornering relationship used in the current version of the 
Von-Schlippe-Dietrich tire model. 


Although Equation (15) does not appear to resemble the Moreland cornering relation- 
ship, it can be shown that they are similar. If the infinite series representing 
the distance delay is restricted to the first two terms, 1 - (h/v)D , and the fol- 
lowing kinematic relationship is introduced: ^ 


= -®t(ya + Rg0B - 


(16) 


where V is the forward velocity of the airplane, then Equation (15) can be 
rewritten as: 


/I + hc>,\ u / . (i + hc^\ . \ 


(17) 


Since the side force is given by 


+ c^A 


(18) 


Equation (17) can be written 


K^(l + hC;^ ) 


Ct • C,. h(Ki + CTDf) 

F„ - +-L — I — lA- 

•'i 


VKj^(l + hC^ 


t) / • 


(1 + Cy^h) 


(19) 


For large values of V the far right hand term is small. Thus Equation (19) is 
similar to Equation (1) with 


( 1 + hC^) 


~ C or — from Equation (5) 
N 


(20) 


and 


( 21 ) 
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From the latter it is apparent that Moreland's cornering relationship is equiv- 
alent to an approximation of the Von Schlippe-Dietrich model. This explains why the 
two models give similar results when used in shimmy analysis. It has been found 
that in many cases the Moreland model is slightly more unstable than the Von 
Schlippe-Dietrich model. It is often more conservative to use the Moreland model 
for landing gear design. Moreland's model is also more simple from a computational 
standpoint since one does not have to deal with the infinite series or, equivalently, 
the time delay function represented by Equations (14) and (13) respectively. 

The forces acting at the tire footprint are only part of the shimmy problem 
formulation. Figure 2 illustrates the factors considered in the total problem, 
namely the airplane, the local structure between the airplane and the landing gear, 
the landing gear structure, and finally the footprint of the tire. 

The fuselage is represented by the three dimensional motion of a fuselage 
reference point. The three motions are a lateral motion yf, a roll motion 4>f, and a 
yaw motion;//^. Three input forces corresponding to these three directions act from 
the landing gear. The motion of the fuselage is given by solution to model equa- 
tions covering the first several normal modes of the fuselage. The equation of 
motion of the fuselage is: 


[m] j*qj + (1 + i 77 f ) [k] jqj = [ 5 ] 


T \]ya 
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where 

[M] = fuselage modal mass matrix (kg) 

[Kl = fuselage stiffness matrix (N/m) 

= loss factor for fuselage structure 

[6] = modal deflection matrix 

forces and moments in the y^xf/and 0 directions 

q = vector of n normal modes 



To obtain the fuselage deflection at the reference point the q vector is multiplied 
by the transpose of the modal deflection matrix. 



The relative deflections due to the structure from the fuselage reference point to 
the landing gear attach point are given by: 
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where 


^t-^f 




r= relative deflections in the y, i// and 0 directions 


[f] = flexibility matrix 


y 

V 

T0 


forces and moments in the y, \}/ , and 0 directions 


The mathematical model of the landing gear is fairly straightforward. Care must 
be taken in properly accounting for cross coupling of flexibilities and inertial 
cross coupling brought on by off-center masses of the landing gear system such as 
torque arms, steering cylinders, lights and attachment fittings. Since each gear 
differs in details the equations of motion will not be given here. 

The remaining equations are those of the tire forces and kinematic restraints. 
For both the Moreland and Von Schlippe-Dietrich models, these equations are: 


F^. = K A . + n A . (25) 

nj 1 j L j 

for the side force, 

^ (26) 

for the moment, and 

V( 0' - = -y^ - Rg0g + Aj (27) 

for the kinematic constraint. 
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For the Moreland tire model the tire equations are completed using Equation (1) as 
the final relationship while for the Von Schlippe-Dietrich model Equation (15) is 
used. 


For dual wheel landing gear systems there is an additional moment acting from 
the ground about an axis normal to the ground plane. This moment, M^^, is governed 
by the following differential equation: 


M + (c R R„/VI ) M + 
w w r g p w 


2K /I /m ( 

w P / w 


dt 


(C„l//2V) - K„C„l//VI (28) 


where in addition to previously defined parameters 

2 

Ip = polar moment of inertia of one wheel and tire (kg-m ) 

Ig = distance between v/heel centerlines (m) 

= torsional spring rate of corotational axle (N-m) 

If = 0 then Equation (28) still may be used by simply setting = 0 and retain- 
ing the remaining terms. It can be noted that an additional tire parameter, is 
needed to deal with dual wheel landing gear systems. This parameter defines the tire 
drag forces in terms of the percent tire slip for small slip ratios. 

Figure 3 shows the Space Shuttle nose landing gear (NLG), which was analyzed 
extensively for shimmy stability using computerized solutions of the equations of 
motion previously described, dynamometer testing, and "runway" tests on the Langley 
Landing Loads Track. This work was extremely important since there were no high- 
speed taxi tests in the Space Shuttle development program. The first "test" of 
the landing gear was an actual landing at over 200 mph. 

The computer study was used to select the optimum damping for the steer damp 
unit, to determine sensitivity to wear, friction, and tire parameters, and generally 
to establish what the margin of stability was for the landing gear system. This was 
then verified and improved upon by the dynamometer test. Application to a flat 
runway was verified by comparison of dynamometer and Langley Landing Loads Track 
test results. The Langley Landing Loads Track could not be used for the full range 
of speeds and vertical loads because of limitations on track capacity. The dyna- 
mometer tests covered the full range of speeds and loads. 

Figure 4 demonstrates, for a typical case, the correlation between the computer 
analysis and the dynamometer (120" flywheel) results. Figure 5 shows the correla- 
tion between the dynamometer tests and the NASA dry concrete runway data from the 
Langley Landing Loads Track. Both sets of results correlate closely, showing that 
the analytical model can be used with confidence for prediction of the Space Shuttle 
NLG's stability over a wide range of landing and roll-out conditions. 

The process of test and analysis for shimmy stability is described in Reference 

7. 
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MAIN LANDING GEAR SHIMMY 


In addition to NLG's, shimmy can be encountered on main landing gear systems 
(MLG's). Three very popular commercial jet aircraft having dual wheel main landing 
gears experienced similar MLG shimmy problems when they were first introduced. The 
problem was traced to a side brace configuration which caused twist in the landing 
gear due to application of a side load. The side braces were all connected to the 
forward part of the shock strut for the three different aircraft. Modifications 
were made to the MLG's to stabilize the systems. In two cases a torque arm damper 
was introduced while in the third case a negative (forward) mechanical trail was 
used to stabilize the shimmy. 

In general MLG's are more complex than NLG's; however, the same tire para- 
meters are needed for MLG stability studies as are used for NLG's. The complexity 
is associated with the landing gear structure itself. This is illustrated in 
Figure 6 which is a drawing of the six wheeled C5A MLG. The analysis of this system 
used a 48 degree of freedom analysis reduced to a twelve degree of freedom system 
for non-linear simulation studies. The 48 degrees of freedom were necessary to 
develop an accurate mass and stiffness matrix for the reduced system. 

Sample stability results are given in Figure 7 and illustrate the importance of 
accurate tire parameters. Here the only difference between the two runs is the yaw 
time constant. The solid curve corresponds to a yaw time constant (C^^) of 0.06 
seconds while the dotted curve applies to a Cj value of 0.03 seconds. The actual 
yaw time constant for the 49 x 17 tires is considerably lower than both of these 
values but the comparison served to show the sensitivity of shimmy stability to tire 
parameters. 



Figure 6.- C5A main landing gear. 
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The net result of the shimmy stability study showed the C5A MLG to be stable, a 
result which was later verified by in-service experience. The analysis was impor- 
tant since this was the first time that this six wheeled landing gear configuration 
was used on an aircraft. 

The diversity of MLG designs is illustrated by the complex single wheel MLG 
shown in Figure 8. This was a preliminary design for the gear which was complicated 
by the fact that it had to retract around a wing mounted rocket. Needless to say no 
single set of equations of motion suffice to cover all landing gear systems. 
Fortunately, however, the tire parameters previously given apply to all landing 
gears. 



'□on .K .10 .IS .» .JO .K .no .15 .iO 

Time, Second 


Figure 7.- Comparison of transient motion for two values 
of the yaw time constant (C5A MLG). 



Figure 8.- Single wheel MLG design for a fighter aircraft. 
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BRAKE WHIRL 


Brake whirl is a highly destructive self excited vibration of the wheel, tire, 
brake and landing gear structure. It is caused by a localized high pressure region 
in the brake stack that brings about out of plane bending of the axle. The result- 
ing motion is characterized by a sinusoidal bending of the axle which takes place in 
two planes, the z-y plane and the x-y plane, as shown in Figure 9. The two motions 
are 90° out of phase which makes it appear that the bending motion is whirling about 
the axle centerline at a speed corresponding to the radial natural frequency. 

The mathematical analysis of the brake whirl problem requires that a number of 
tire properties be known. Since the frequency of the whirl motion is generally in 
the range of 200 to 400 Hz for most landing gears the tire does not respond as it 
does for low frequency motions. The tire tread remains fixed in the x-z plane except 
for its steady forward motion and rotation. The heavy bead area of the tire moves 
with the wheel adding to its effective mass and its moments of inertia. The region 
between the tread and the wheel is treated as a spring with radial and angular 
spring constants and associated loss factors. 


y 


Figure 9.- Brake whirl. 
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The specific tire parameters are: 


K = spring rate of tire associated with radial displacement of the wheel 
with the tread held fixed (N/m) 

7]^ =. loss factor associated with 

= diametral angular spring rate of the tire associated with twist of the 
wheel with the tread held fixed (N*m/rad) 

\r ~ factor associated with 

Mtw = tnass of the tire moving with the wheel (the bead area of the tire) (kg) 
Ip = polar moment of inertia of m^^ (kg-m ) 

'2 

Idt = diametral moment of inertia of (kg*m ) 

= tire rolling radius (m) 


Brake whirl stability is a somewhat more serious problem on carbon brake systems 
because of their mass distribution properties. A careful analysis as well as design 
changes to stabilize the motion are important in the initial stages of the brake 
development. 

Measurement of the tire stiffness parameters is difficult. As a first approx- 
imation is taken to be the average of the fore and aft tire stiffness and the 
vertical stiffness of the tire. The value of is calculated based on lateral de- 
formation of the tire at both the footprint and 180° from the footprint. The loss 
factor is estimated to be approximately 0.15 for typical aircraft tires. The mass 
Mtw ''s roughly 30% of the tire for typical aircraft tires. Calculation of Ip and 
I(jt can be carried out knowing M-^w and the location of the bead bundle. The^rolling 
radius can be determined by conventional methods and is equal to the undeflected 
radius of the tire minus one-third of the tire deflection as a very close approxi- 
mation. 
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BRAKE CHATTER AND SQUEAL 


Brake chatter and squeal are self excited motions of the brake and landing gear 
system brought on by so called negative damping at the rotor-stator interface. The 
negative damping results from an increase in brake torque for a decrease in the in- 
stantaneous brake slip velocity. This characteristic yields a term in the equations 
of motion which looks like a damping term in that it is a force or torque propor- 
tional to velocity, but has a sign opposite to that of the dissipative damping of the 
motion. Thus, the lining characteristic has been termed negative damping. 

The chatter motion is a low frequency motion (5 Hz to 15 Hz) and consists of 
fore and aft bending motion of the shock strut. Squeal consists of wind-up of the 
stationary parts against the torque take-out system and has a natural frequency of 
150 to 500 Hz. 


A relatively simple model of brake chatter and squeal uses three degrees of 
freedom. These are the fore and aft bending motion of the landing gear shock strut, 
the torsional wind-up of the nonrotating parts of the brake against the torque take- 
out system, and the forward velocity of the airplane. The nonrotating parts of the 
brake consist of the brake piston housing, the pressure plate, the backing plate, 
the stators of the brake, and the torque tube of the brake which carries the 
stators. The torque take-out system consist of an axle flange on most single and 
dual wheel landing gears and a brake equalizer rod on truck type landing gears com- 
monly used on larger aircraft. 


The tire parameters used in the simplified chatter and squeal analysis consist 
of the following: 

R^ = tire rolling radius (m) 

Rg = loaded radius of tire (m) 

= slope of drag force versus slip ratio for small slip ratios (N) 

= mass of tire (kg) 

2 

Ip = polar moment of inertia of tire (kg-m ) 

C = ratio between the shift in the center of pressure and the aft deflection 
^ of the footprint 


This three degree of freedom analysis is quite simplified but is useful for 
preliminary studies of stability of chatter and squeal, particularly for the 
selection of lining materials that will not result in self-excitation of the vi- 
bratory motions. At a later stage in the landing gear development, a more detailed 
analysis of brake chatter and squeal such as the one that is described in the next 
section can be carried out. 
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COMPLEX COMBINATIONS OF CHATTER AND SHIMMY 


For most modern main landing gear systems the chatter motion is not simply a 
fore and aft bending action of the shock strut but rather consists of a complex 
motion of the gear involving fore and aft, lateral, roll, and yaw motions of the gear. 
Thus many of the same tire parameters that are used in shimmy analysis are also 
needed for an analysis of the brake chatter and squeal. 

Figure 10 shows a typical truck type main landing gear system and some of the 
nomenclature and motions considered. There are ten degrees of freedom in the 
model : 

Xg = aft deflection at the truck pivot point (m) 

Vg = lateral deflection at the truck pivot point (m) 

0g = roll angle of truck (rad) 

0 g = end bending angle of shock strut (rad) 

0g = yaw angle of truck (rad) 

= angular velocity of one of the forward wheels (rad/sec) 
cOg = angular velocity of one of the aft wheels (rad/sec) 

^sl ~ angular rotation of the forward brake stationary parts (rad) 

©ss = angular rotation of the aft brake stationary parts (rad) 

Xp = forward translational motion of airplane (m) 

A fairly extensive set of tire parameters is needed for the model. The shimmy 
action uses the simple cornering relationship equating the tire lateral force to a 
sideslip angle divided by the coefficient of yaw. This is equivalent to the modified 
Moreland model with C^ = 0. 
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The tire parameters used in the model are: 

= fore and aft tire spring rate (N/m) 

= vertical tire spring rate CN/m) 

C = coefficient of yaw (rad/N) 

lifP = lateral coefficient of friction between tire and runway, a function 
of slip ratio 

^\f/m ~ niaximum value of 

11 ^ = ratio between self-aligning torque and sideslip angle (N-m/rad) 

= maximum fore and aft coefficient of friction between tire and runway 

C = ratio between the drag force and the slip ratio for small slip 
ratios (N) 

= ratio between the shift in the center of pressure and the aft 
deflection of the tire 

= mass of tire (kg) 

2 

= polar moment of inertia of inflated tire (kg>m ) 

2 

Idt = diametral moment of inertia (kg-m ) 

Rg = loaded radius of tire (m) 

= rolling radius of tire (m) 


A sample result of the nonlinear analysis is shown in Figure 11. In this run 
the negative damping from the brake lining was set relatively high and the system 
damping, made up of several nonlinear terms, was set at about one half the expected 
value. A pulse was added to the torque to cause excitation at the midpoint of the 
10 mph stop. It can be seen in Figure 11 that the initial motion is stable but 
that the pulse causes enough disturbance that the resulting motion is near neutral 
stability. The mode shape of the motion is far from a simple fore and aft truck 
motion. Note that the lateral motion of the truck is the same order of magnitude as 
the aft motion of the truck. Also, there is some yaw motion and considerable roll 
motion of the truck. 

Studies of this type help to establish the amount of negative damping that can 
be tolerated by both the chatter motion and the high frequency squeal motion of the 
landing gear. Also mode shapes and frequencies of the several low frequency chatter 
motions can be established together with the margin of stability of each mode. 
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Figure 10.- Nomenclature and motions for complex chatter 
and shimmy analytical model. 
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Figure 11.- Sample computer run for MLG chatter and shimmy analytical model 



OTHER ANALYSIS AREAS 


The other analysis areas mentioned at the beginning of this paper also require 
detailed information on tire characteristics. Anti-skid performance requires a 
model of the brake pressure to brake torque transfer function along with the 
dynamic characteristics of the anti -skid system. The remainder of the model is 
similar to that used for chatter and squeal. Gear walk which is a stick-slip 
motion similar to chatter can come about due to stick-slip action at either the 
ground surface or the brake interfaces. Rough terrain operations require a tire 
model capable of predicting vertical and drag loads due to obstacle or hole impact 
at obstacle lengths which range from a small fraction of the footprint length up to 
several times the footprint length. The short lengths present more of a problem but 
methods are available for predictions of tire loads for all lengths. Rough terrain 
operations analysis also requires a detailed model of the shock strut stroking 
dynamics in addition to the tire model. 


CONCLUSIONS 


Tire dynamic parameters play a key part in many design analysis areas for air- 
craft landing gear systems. Advanced analytical methods to predict tire parameters, 
such as finite element methods, are needed together with experimental methods for 
verification of predictions for selected cases and determination of parameters that 
cannot be determined analytically. Accurate tire data allows more accurate design 
analyses and correct prediction of dynamic performance of aircraft landing gear. 

The net result is a more reliable and more efficient system. 
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